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This  is  a  working  paper.  It  may  be  expanded, 
modified,  or  withdrawn  at  any  time.  The  views, 
conclusions,  and  recommendations  expressed  herein 
do  not  necessarily  reflect  the  official  views 
or  policies  of  the  United  States  Air  Force. 


Technical  problems  eisaociated  vlth  hardened  base  facilities  have  been 
under  continual  Investigation  at  RAND.  A  part  of  this  work  concerns  the 
design  of  shelter  structures  to  withstand  the  blast  or  pressure-rise 
effects  of  a  nucletfr  explosion.  This  research  memorandum  treats  one 
aspect  of  the  problem. 

A  generalized  mathematical  model  Is  presented  of  the  rise  of  pressure 
with  time  within  an  enclosure  as  a  function  of  enclosed  volume,  leakage¬ 
opening  area,  static  overpressure,  and  nuclear  weapon  yield.  The  results 
are  shown  graphically.  Their  applicability  to  the  design  of  protected 
enclosures  is  illustrated  by  two  exan5)lea  of  personnel  shelters:  a  nmAi i 
home  shelter,  and  a  considerably  larger  arremgement  housing  up  to  a 
hundred  people . 

This  Information  should  be  useful  in  evaluating  protective -construction 
design,  particularly  as  It  is  affected  by  the  incorporation  of  heating 
and  ventilation  subsystems . 
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SUMMARY 

A  blast  shelter  or  other  sliBllar  enclosure  normally  has  numerous 
openings  leading  to  the  outside  for  Ingress  and  egress,  ventllatloQ,  and 
utility  lines.  Thus,  there  Is  a  possibility  of  air  leakage  because  of  an 
Ineffective  or  defective  door  seal,  or  an  Inoperative  or  open  valve  In  the 
ventilation  system. 

Any  leakage  to  the  Inside  of  high-pressure  air  resulting  from  a 
nuclear  blast  wave  would  tend  to  raise  the  pressure  within  the  shelter. 

The  magnitude  of  the  rise  depends  on  the  shock  strength,  the  time  duration 
of  outside  pressure,  the  sizes  of  the  leakage  area,  and  the  shelter  volume. 

For  certain  conditions,  such  as  large  Internal  shelter  volume  and 
small  opening  areas,  such  Items  as  door  seals  and  blast  valves  In  the 
ventilation  system  may  not  be  necessary.  Under  these  conditions,  a  nuclear 
blast  wave  could  not  cause  an  appreciable  pressure  rise  within  the  enelosiire 
without  physical  destruction  of  some  part  of  the  structure. 

A  somewhat  Idealized  model  of  a  blast  shelter  has  been  analyzed  for  a 
large  range  of  the  applicable  variables.  The  pressure-time  history  and  the 
peak  value  of  the  Internal  pressure  have  been  determined.  The  results  of 
the  analysis  are  presented  In  generalized  form  by  means  of  a  series  of 
curves.  In  addition,  two  examples  are  presented  which  are  representative 
of  personnel  blast  shelters:  first,  a  typical  home  shelter,  and  second, 
a  shelter  that  might  serve  50  to  100  people. 
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SYMBOLS 

A  effective  orifice  area  or  leakage  area 

B  constant 

C  velocity  of  sound 

D  constant 

E  internal  energy 

g  gravitational  constant 

H  enthalpy,  y  E 

J  mechanical  equivalent  of  heat 

K  constant  ■  y 

P  absolute  pressure 

q  dynamic  pressure 

R  perfect  gas  constant 

^  absolute  air  temperature  in  degrees  Bankine 

T  absolute  temperature 

t  time 

t^  time  duration  of  positive  phase  of  overpressure 
V  shelter  volume 

W  weight  of  air 

W  weight  flow  rate 

X  pressure  ratio  p2./^21 

pressure  ratio 
y  pressure  ratio  P2/*’21 

CX  non-dlmunsloual  time  t/t 

o 
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^  constant  In  exponent 

Y*  ratio  of  specific  heats  «  l.i»  for  air 

A  Increment 

AP^  Initial  overpressure  •  -  1 

&  constant  In  exponent 

0  sonic  fill  time 

V  velocity 

p  weight  density 

00  function  of  ( ) 

Subscripts 

1  outside  conditions 

2  Inside  conditions 

1  Initial  conditions 

0  peak  Initial  conditions 

B 


aaxlinuB  value 
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I.  IWHODUCTIOW 

In  the  design  of  blast  proteotiTa  shelters  or  other  enelosures,  it  is 
iBg)ortant  that  attention  he  giren  to  door  seals  and  ventilating  openings. 

One  reason  is  that  the  structural  loading  of  an  enclos\ire  nay  be  a  function 
of  the  pressure  rise  resulting  froa  leakage,  as  in  the  case  of  open  hatches 
on  a  ship's  hold. 

An  enclosure,  when  subjected  to  an  exterior  blast  vare  of  finite 
duration,  will  leak  froa  any  opening.  The  pressure  vlthin  will  rise  to 
soae  peak  value  and  then  decay  as  the  flow  through  the  opening  reverses. 

Door  seals  and  check  valves  in  ventilating  systeas  aay  be  necessary  to 
Halt  the  Inside  pressure  to  a  tolerable  levelj  coirrespondlngly,  an  open 
hatch  alloying  pressure  equalization  aay  protect  the  deck  or  hull  of  a 
ship  froa  failure. 

The  purpose  of  this  stvuly  is  to  analyze  the  pressure-tlae  history 
vithln  an  enclosure  as  a  function  of  enclosed  volune,  leakage  opening 
area,  static  overpressure,  and  nuclear  veapon  yield.  The  peak  pressure 
vithln  the  shelter  for  any  specific  set  of  conditions  aay  then  be  deteralned. 
This  Information  should  be  useful  In  engineering  design  as  veil  as  In  the 
analysis  of  blast  resistance. 
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II.  ASSUMPTIOWS 

Tb«  Mclosed  voluma  or  shelter  Is  e  single  cosqmrtaent  having  a 
volune,  Vgi  and  an  opening  or  orifice  of  area,  A,  leading  to  the  outside. 
Further,  the  area.  A,  is  siiall  cosipared  to  the  voluise}  that  is,  the  area 
■ay  be  Measured  in  fractions  of  a  square  foot,  while  the  volume  may  be 
described  by  hundreds  of  cubic  feet. 

nie  shelter  volume  opening  system  is  critically  damped;  that  is,  the 
cavity  does  not  resonate  as  a  result  of  sin  external  pressure  stimulus. 
Viscous  losses  from  ojiening  auod  furnishings,  other  losses,  and  the  limiting 
effect  of  the  normal  shock  wave  choking  the  opening  and  regulating  the  flow 
during  the  period  when  most  of  the  dabs  flow  occurs  tend  to  inhibit  oscil¬ 
latory  pressure  fluctuations.* 

The  initial  absolute  pressure,  within  the  volume  or  enclosure 

is  ainbient  atmospheric  pressure  (corresponding  to  IU.7  pel  at  sea  level) 
prior  to  the  arrival  of  the  air  shock  wave.  Upon  arrival  of  the  air  shock 
the  exterior  portion  of  the  oi>ening  experiences  a  sharp  pressure  rise. 

In  this  analysis  the  peak  exterior  absolute  pressure,  ^1'  is  assumed  to 
be  the  sum  of  the  ambient  pressure,  Pg^,  and  the  pecUc  static  overpressure, 
Ap,  at  the  shock  front.  This  msy  or  may  not  be  the  proper  estimate  of 
the  external  pressure,  depending  on  the  geometry  of  the  shock  front  and 
the  exterior  surface  of  the  enclosure  containing  the  opening.  Tbus  in 
any  calctilation  it  nay  be  necessary  to  modify  the  value  of  the  overpressure, 
AP.  As  an  example,  consider  the  case  of  a  door  jamb  leakage  where  the 
door  is  flush  with  the  exterior  of  the  enclosure  which  is  parallel  to  the 

piano  of  the  arriving  air  shock  wave;  that  is,  the  shock  strikes  head-on 
# 

Despite  these  damping  effects,  the  criterion  of  resonance  should  be 
considered  in  shelter  design. 
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with  the  opening.  In  thie  case  the  peak  forcing  preeaurei  would  be 
the  sum  of  the  preshock  ambient  pressure  the  peak  static  overpressure 
AP|  the  diffraction  pressure,  and  the  peak  dyneunle  pressure. (Dlf* 
fraction  effects  are  of  primary  Importance  In  cases  of  very  high  over¬ 
pressures  of  about  100  atmospheres.  The  effect  Is  usually  of  minor 
importance  at  low  pressures  of  a  few  atmospheres,  at  least  as  far  as 
leakage  Is  concerned.)  On  the  other  hand.  If  the  opening  were  on  the 
downstream  side  of  the  structure  the  peak  forcing  pressure  would  be  the 

sum  of  and  /iP  minus  a  constant  times  the  peak  dynamic  pressure 

(4) 

The  value  of  the  constant  depends  on  the  exterior  shape  of  the  body 

(shelter)  and  the  particle  Mach  number,  and  can  vary  from  a  maximum  value 

of  about  1.4  at  low  subsonic  speeds  to  a  value  as  small  as  0.1  at  Mach 

At  tero  euigle  of  incidence,  surface  normal  to  the  plane  of  the  shock  wave, 

the  constant  becomes  zero.  Illustrative  calculations  in  this  memorandum 

will  be  on  this  basis .  This  condition  can  correspond  to  an  opening  through 

the  flat  horizontal  top  of  m  enclosure.  The  outside  pressure  variation 

* 

with  time  Is  described  by  the  relationship 

?!  -  Pgi  MPii  -Pgi)  (1  -  1“)  j^B  e  ■  P^/^o  +  D  e'  (l) 


where  the  peetk  static  ovexTpressure  4P  -  -  p^^  and  B,  ,  D,  suid  A 

depend  on  the  peak  static  overpressure.  The  time  t^,  the  duration  of  the 
positive  pressure  pulse,  is  a  function  of  both  peak  pressure  and  weapon 
yield.  Equation  (l)  can  be  rewritten  in  non-dimensional  form  as 


'b. 

’“21  J 

>  1 

(2) 


See  the  appendix. 
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The  Internal  or  shelter  preeaure,  Pgf  la  aeeumed  to  be  equal  to  the 
static  preasure  of  the  air  flowing  from  the  nozzle.  That  Is^  there  Is  no 
recovery  of  the  dynamic  pressure  In  the  air  stream  from  the  opening. 

Air  acts  aa  an  Ideal  diatomic  gas  In  the  overpressure  range  of 
interest  (zero  to  perhaps  20  atmospheres). 

The  system  is  adiabatic.  It  is  assumed  that  no  energy  is  lost  to  the 
walls  as  the  air  is  forced  through  the  opening.  In  any  real  system  there 
%rould,  of  course,  be  heat  loss,  especially  as  the  effective  length-to- 
diameter  ratio  of  the  opening  increases.  This  assumption  results  in  the 
prediction  of  slightly  higher  internal  pressures  than  would  actually  be 
experienced. 

The  area  of  the  opening  la  assumed  to  be  the  mlnlaua  section  area  in 
any  conduit,  seal,  or  other  opening.  However,  in  any  real  ease  there  are 
losses  such  as  the  diffraction  effects  at  the  entrance,  expansions  and 
contractions  within  the  channel,  viscous  effects  at  the  walls,  inertia  or 
vena  contracta  effects  at  the  minimum  section,  and  exit  effects.  These 
effects  all  tend  to  restrict  the  flow;  therefore,  the  effective  area  of 
any  actual  opening  would,  in  general,  be  somewhat  smaller  than  the  mlniaua 
geometric  area.  This  assumption  also  results  in  the  prediction  of  internal 
pressures  somewhat  higher  than  in  any  real  case. 

The  system  is  choked  (flow  velocity  equal  to  Mach  l)  at  the  mlnlmuffl 
section  of  the  opening  as  long  as  the  pressure  ratio  remains  critical 
(^2  <0.528  for  air). 
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III.  AHALYSIS 

Hie  velocities  and  pressures  involved  can  be  described  by  an  energy 
balance  across  the  nozzle.  For  an  ideal  ges  under  adiabatic  conditions 
and  on  a  unit  weight  basis,  the  sum  of  the  enthalpies  and  kinetic  energies 
should  be  equal: 

V  ^  V  ^ 

^  ®1  fg”  ’  J  r  ^2  *  ?|"  ^3) 

However,  the  Initial  velocity  la  very  low  and  nay  be  neglected. 
Solving  for  the  nozzle  exit  velocity  we  have 

=  2gJ  r  -  Eg)  (U) 

The  weight-flow  rate  through  the  nozzle  may  be  written  as 

W  -  Pg  A  Vg  «  pg  A  Wg  J  r  [e^  -  Eg)  (5) 

The  unit  internal  energy  of  a  perfect  gas  can  be  expressed  in 
terms  of  pressure  and  density  as 

The  weight-flow  rate  (Eq.  (3))  may  then  be  rewritten  as 


With  the  assumptions  of  adiabatic  flow,  no  molecular  changes,  and  re- 
^2  /P2\ ^ 

versible  flow,  - —  ,  we  can  describe  the  weight-flow  rats  in  terns 

U  \Pl/ 

of  the  inlet  conditions  and  the  pressure  ratio  across  the  nozzle, 
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Thus  Bq.  (7)  bscooiss 

W  -  A  2g  ^ 


•e 

to 

to 

fPji 

lii-ll 

r 

t - 

1 

1/2 


(8) 


The  flow  z«ts  is  thus  dasoribsd  as  a  function  of  the  Initial  state  of 
the  ehocked  gae  exterior  to  the  ehelter,  the  orifice  area,  and  the  pres¬ 
sure  ratio  across  the  orifice. 

It  should  be  noted  that  the  velght-flov  rate  W  ae  described  by 
Eq.  (6)  passes  through  a  naxlmua  \dieD  plotted  as  a  function  of  the  pres¬ 
sure  ratio  ^2/^1*  pressure  ratio  at  v^lch  this  naxlmia  occurs  Is 
a  function  of  the  ratio  of  the  specific  heats  of  the  gas.  The  isaxisam 
for  any  two  of  the  three  equatlon-of -state  variables  held  constant  can 
be  deteralned  by  setting  d  w/d  Pg,  derived  froa  Eq.  (8),  equal  to  sero 

and  solving  for  the  nozzle  pressure  ratio: 

r 

^  (9) 

For  normal  air  (r  «  1>^)  the  value  of  the  pressure  ratio  at  the  siaxlioua 
flow  rate  is  very  close  to  O.528. 

Recalling  the  energy  balance  froo  Eq.  (U),  the  definition  of  Inter¬ 
nal  energy  (Eq.  (6)),  and  the  adiabatic  relationship  between  density 
and  pressure,  we  have  by  mutual  substitution  an  expression  for  velocity 
in  terms  of  the  nozzle  state  variable  and  the  pressure  ratio: 


M 

2 

IKC 

r?r  if 


in 


-  1 


(10) 
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We  may  now  aolve  for  the  velocity  of  maximum  flow  by  substitution  of 


—  from  Eq.  (9)  in  Eq.  (lO): 

^1/  max 


(11) 


Thus,  the  flow  velocity  at  this  critical  condition  is  equal  to  the 
local  velocity  of  sound  at  the  nozzle  exit,  C^. 

The  nozzle  would  then  limit  the  flow  rate  W  (for  any  specific  set 
of  inflow  state  conditions  and  p^)  to  the  constant  choked  or  sonic 
flow  as  long  as  the  ratio  of  pressure  within  the  shelter  to  the  outside 
pressure  is  less  than  0.528.  The  critical  flow  rate  during  choke  flow 


would  be  represented  by 


^2  ''2m 


r  ^1 

=  ^)2v'2g^  — 


For  pressure  ratios  greater  than  O.528  the  flow  velocity  would  be 

less  than  sonic  and  would  be  described  by  Eq.  (8).  The  solid  portion  of 
the  curve  of  Fig  1  is  descriptive  of  the  instantaneous  flow  rate  cor¬ 


responding  to  any  value  of  the  outside  conditions  and  any  pressxire  ratio 

across  the  nozzle.  The  flow  rate  variable  W  is  that  value  of  air  flow 

c 

described  by  any  two  of  the  three  state  variables  and  T,  and 

values  of  P^  and  P^^  such  that  the  ratio  Is  O.528.  Ibe  ordinate  of 

Fig.  1  is  thus  expressed  in  multiples  of  the  critical  flow  rate.  The  ex¬ 
pression  for  the  curve  between  O.528  and  1  is 
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llie  pressure  and  temperature  within  the  enclosure  may  be  related  by 
the  perfect  ^s  lav 


A  heat  balance  relates  temperature  and  flow 


(lU) 


i”l2  *  ^21  ^2  -  ^  C, 


(15) 


Eliminating  between  the  above  equations  and  solving  for  ?2  ve  have 


r  « 

?2  J'  *  Pgi 


(16) 


The  variation  of  the  internal  pressure  with  time  nay  be  determined  by 
differentiating  with  respect  to  time 


r  R  T 


11 


12 


- 

\  ^c/ 


(17) 


'2  '2 

The  ratio  w/w^  and  the  expression  for  have  been  detenained  previously. 
Substituting  the  expression  for  W  from  Eq.  (12)  in  Bq.  (17)  we  have 


^  “^li  .  !  „  _  /  2 


^2  “  ^  i^lPl®  nlj 


r-1 


w 


12 


(18) 


Hi 


But  recalling  that  ~  ~  =  I  pi-j  ^  and  that  -  VfB  RT 


rAOiiPj^  /■’ll 

^  hj  ^ 


lil 

2  r-1 

r+i 


w 


12 


w„ 


11 

(19) 

(20) 
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or  in  non-dlmenaional  fom 


and  0  -  Vg/A 

The  quantity  x  la  defined  by  awana  of  Bq.  (2)  in  dliaenaianleia 

fom 

r  -pa  -  aa, 

X  -  1  +  (x^  -  1)  (1  -  a)  [b  e  4  D  e  j  (22) 

We  now  vleh  to  determine  y  aa  a  function  of  a  (ahelter  preaaure  as 
a  function  of  time)  for  various  values  of  x^  (static  overpressure),  for 
a  rsinge  of  the  variable  (the  ratio  of  the  poaltivo  duration  of  the 

overpressure  to  the  minimum  possible  fill  tine),  and  for  values  of  the 
constants  B,  3 ,  D,  and  6,  which  are  functions  of  the  peak  static  over¬ 
pressure  . 

The  solution  is  valid  only  for  flow  into  the  shelter  or  until  the 
pressure  Inside  rises  to  its  maximum,  which  is  also  equal  to  the  instan¬ 
taneous  outside  pressure  i.e.,  -  Pg-  This  condition  also  corres¬ 

ponds  to  both  W  and  P^  equal  to  zero.  Therefore,  as  shown  by  Eq.  (21) 
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and  Fig.  1,  for  this  coodltlon 


sad 

r  -aa  -aai 

-  X  -  1  +  (x^  -  1)  (1  -  a)  [Be  ♦  D  e  J  (24) 


Thus  we  knov  that  the  loci  of  oaxina  In  y  for  any  discrete  value  of  x 
■ust  fall  along  the  curve  of  Eq.  (24).  This  la  illustrated  In  Fig.  2. 


F'9.2 —  Loci  of  solutions  for  maximum  shelter  pressure 

os  0  function  of  time 
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Of  course  no  solution  for  the  naxlisum  exists  at  auay  tins  greater  than  t^, 
the  duration  of  the  posltlYs  phase  of  the  overpressure. 

This  does  not  provide  an  explicit  solution;  in  fact  there  Is  an 
infinite  nunber  of  solutions  for  the  naxljna  Blaag  the  curve  of  Fig.  2. 

The  problea  of  Integrating  Eq.  (21)  reaains. 
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IV.  REoULTS 


Equation  21  has  be«n  solved  numerically  for  several  values  of  the 
static  overpressure:  1,  5#  and  10  atmospheres.  The  results  are 

presented  In  the  curves  of  Figs.  3  to  6.  The  calculations  of  internal 
pressure  as  a  function  of  time  have  been  carried  only  to  the  condition  of 
maximum  internal  pressure.  In  general,  the  pi'cssure  decay  times  will  be 
greater  than  the  rise  times. 

Figure  7  Ifl  a  summary  of  the  generalized  solution  for  the  maximum 
pressure  rise  to  be  expected  within  the  shelter. 
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Fig.  6 — Pressure-time  history  within  on  enclosure 
Xq»  II,  ■  10  otmospheres 
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V.  NIIMKHICAL  EXAMl’LE.i 


Example  1 

As  an  extunple,  let  us  consider  a  shelter,  representative  of  a  small 
home  blast-ond-fiillout  shelter,  with  the  following  significant  proiiertles: 


Volume 

Effective  leakage  area 

(3 -In.  Inlet  and  3 -in.  exliuust  openings) 


Vg  =  500  ft^ 
A  =  0.1  ft^ 


WeapCn'yleld 


3-t'fr  surface  burst 


Peak  static  overpressure 


aP  »  3  atmospheres 


Since  the  Initial  pressure  ratio. 


21 


11 


jj-  ,  is  suiorcritical  (less 


than  0.53),  at  the  onset  of  the  shock  front  we  are  Interested  In  the 
temperature  and  the  sonic  velocity.  The  peak  temperature  associated  with 
the  shock  front  may  be  read  from  Fig.  9  (p*  29).  The  value  for  AP  =  3 
atmospheres  Is  830°R.  The  Initial  choked  velocity  through  the  orifice  Is 

then 


C 


11 


1120 


820 


520 


1415  ft/sec 


where  the  assumed  preshock  ambient  temj^irature  and  sonic  velocity  are  520  U 

and  1120  ft/sec,  respectively. 

The  sonic  fill  time  9  then  becomes 
V„ 


9  = 


AC 


11 


500  ^ 

o.r(iiii5)  ° 


Erode *3  values  of  the  positive  duration  time  of  the  static  overpressuii 
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pvilse  are  presented  In  Fig.  8  (p.  27)*  The  values  are  descriptive  of 
a  Bea~level  1-KT  nuclear  explosion  which  has  been  alrburst.  The  positive 
duration  of  the  alrburst  at  an  overpressure  of  3  atmospheres  Is  seen  to  be 
0.093  see.  Assuming  perfect  reflection  for  a  surface  burst  of  a  3-MT 
weapon,  the  duration  Is 


tg  -  .093  (2)^/^  (3000)^/^  «  1.73  sec  (3  MT) 
Thus  the  time  parameter  becomes 


t 

0 

0 


hJl 

3.53 


-  0,49 


The  solution  for  the  maximum  Internal  pressure  may  be  read  from 
Fig.  7«  Entering  with  the  static  overpressure,  3  atmospheres,  to  a  vaLlue 
of  the  time  parameter,  0,49,  one  reads  a  value  of  peak  pressure  ratio 


«  1.6 


Thus  the  peak  pressure  rise  within  the  shelter  Is  0.6  atmosphere  or 
about  9  psi.  Actually  the  pressure  rise  would  probably  be  somewhat  less. 


Example  2 

If  the  shelter  were  larger,  say  5000  ft^,  had  the  same  leakage  area# 

2 

0.1  ft  ,  and  were  exposed  to  an  overpressure  of  10  atmospheres  from  a  10-MT 
surface -burst  nuclear  weapon  we  would  have: 

Vg  =  5000  ft^ 

A  -  0.1  ft^ 

W  =  10  MT 


AT  =  10  atmospheres 
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The  peak  static  shock  temperature  would  then  be  about  10^0^  and  the 
shock  velocity  woxUd  be 


-  1120 


1590  ft/sec 


The  sonic  fill-time  becomes 


®  ”  o'i^96T  “ 


The  positive  duration  time,  from  Fig.  9t  would  be  0.06  sec  for  a  1-KP 
airbxirst  weapon.  The  duration  time  for  the  surface -burst  10-MT  weapon 
woxild  be 


t^  -  0.08  (2)  (10,000)^/^  -  2.2  sec 


The  time  parameter  would  be 


0"  -  IT  '  *0^1 


The  maximum  pressvure  ratio  within  the  shelter  would  be  (from  Fig.  7) 


2m 

2i 


1.25 


The  peak  pressure  within  the  shelter  would  be  about  3*7  pei.  As  one  would 
suspect,  the  sonic-choke  principle  of  valving  becomes  most  effective  for 
large  shelters  having  long  sonic  fill  times  as  described  by  small  values  of 
the  time  x^arometer. 

There  are  some  complications  evident  in  the  second  example .  In  a 


j 
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ventilating  ayetem  having  a  short  pipe  between  the  exterior  and  the  interior 
of  the  shelter,  hot  luminous  gas  may  pour  from  the  opening.  This  possibility 
can,  of  course,  be  obviated  by  thoughtful  design  of  the  system.  A  good 
description  of  the  environment  at  a  close-in  distance  corresponding  to 
high  static  overpressure  from  a  nuclear  weapon  may  be  obtained  from 


Refs.  1-3. 


Appendix 

PROPERTIES  AT  AND  FOLLWINC  THE  SHOCK  FROffT 
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STATIC  OVERPRESSURE 

The  variation  of  static  overpressure  with  time  depends  on  the  dlstaince 
from  the  explosion  center,  hence  the  i»eak  static  overpressure.  The  over¬ 
pressure  profiles  as  well  as  the  peak  Instantaneous  pressure  have  been 

(l  2) 

described  by  Erode  '  by  means  of  a  dichotomy  of  exionentlals  representing 
two  decay  rates 


where 

+  0.23 


2.28  8  +  AP^^ 


27.66  +  +1.2  +  .007 


D  =.  1  -  B 


-I 

V  1  +  .1 


11 


dp 


li 


d  P 

1.5  11 


1500  +  AP 


3/2 

11 


5  =.  9.0  +  1.4  AP 


and  for  1  KT 


0.1  AP 


11 


°  70-2  AP^^  +  AP^^ 


^  20  AP^ 


(For  >  13  omit  the  natui-ul  logarithm  term.) 
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For  weapon  yields  Y  other  than  1  ICT 

to  =  ^o-lKT  ^  kilotons) 


Values  of  B,  0  ,  T),  S  ,  artd  to  are  presented  In  Fig.  8. 


DYNAMIC  rRE.ioURE 

Although  values  of  the  dynamic  pressure  have  not  been  used  in  the 
calculations  presented  in  this  memorandum,  Brode's  values  are  presented  for 
the  sake  of  completeness.  It  is,  of  course,  desirable  to  include  the 
effect  of  dynamic  pressure  in  the  forcing  pressure  signal  1'^  when  the 
geometry  and  the  consequent  pressure  recovery  demand  it. 

Expressions  similar  to  those  representing  static  pressure  have  been 
used  to  represent  the  variation  of  dynamic  pressure  with  time,  although 
the  dynamic  pressure  decays  much  more  rapidly  and  goes  through  zero  in  a 
shorter  period  of  time,  t^. 


.oq 


11 


E  e 


«t/t^  -  st/t 

+  E  e 


where 


t 

X 


.0421  In 


=  .04 

ip  ■ 

10,000  11 

10,000  +  Ar?, 


(for  <  33) 

(for  >  33) 

(for  >  1) 


E=  1  -  E 


fh>2660 

27 


Fig.8 — Values  of  parameters  describing  peak  pressure  vs  time 
from  a  scaled  l-KT  sea- level  nuclear  detonation 
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«  -  0.25  +  3 .6  V^Pii 
f-T*8V5^  ♦ 


AIR  TEMPERATURE 

The  velocity  through  the  orifice  during  the  period  of  choked  flow  Is 
determined  by  the  local  sonic  velocity,  hence  by  the  local  temperature. 
Hie  functioned  relationship  for  the  sonic  velocity  Is 


'li 


1120 


ft/sec 


where  Is  measured  In  °R  and  assuming  a  60°F  or  520°R  day. 

In  regions  of  low  overpressure,  ^  atmospheres,  the  peak 

temperature  occ\irs  at  the  shock  front  and  may  be  estimated  from  the 
adiabatic  relationship 

^  for  0  <  AP  <  10  atmospheres 


where 


-  520°R  for  a  normal  60°F  day. 

In  regions  of  greater  overpressures  It  Is  desirable  to  compute  tempera¬ 
tures  on  the  basis  of  a  normal  shock.  Values  of  the  peak  static  temperature, 
based  on  a  60^F  or  520^  day,  have  been  computed  and  are  presented  in  Fig.  9* 


Overpressure  (otmospheres ) 


Fig. 9 —  Peak  shock  front  static  temperoture 

(Preshock  ambient  temperature  of  520®R,60®F) 
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